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Financial time series exhibit two different type of non linear correlations: (i) volatility autocorrela- 
tions that have a very long range memory, on the order of years, and (ii) asymmetric return-volatility 
(or 'leverage') correlations that are much shorter ranged. Different stochastic volatility models have 
been proposed in the past to account for both these correlations. However, in these models, the 
decay of the correlations is exponential, with a single time scale for both the volatility and the lever- 
age correlations, at variance with observations. We extend the linear Ornstein-Uhlenbeck stochastic 
volatility model by assuming that the mean reverting level is itself random. We find that the re- 
sulting three-dimensional diffusion process can account for different correlation time scales. We 
show that the results are in good agreement with a century of the Dow Jones index daily returns 
(1900-2000), with the exception of crash days. 
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I. INTRODUCTION 



The interest for financial markets models has considerably increased during the last twenty years. The traditional 
Black-Scholes Gaussian model is now well known to be inadequate; devising faithful and convenient models for the 
price processes is essential for risk control, derivative pricing and other trading operations. Data analysis has now 
become a mandatory step in model building, using the wealth of available tick-by-tick data or high-frequency data . 
There is now enough data to establish, with statistical accuracy, some fundamental properties of market dynamics. 
Empirical observations provide a list of universal features that seem to hold in most markets, quite independently of 
the observation period and of the type of asset. These are the so-called stylized facts (see e.g. 0,0,0), which should 
be reproduced, at least to some extent, by any reasonable market model. 

One of the best established stylized fact is the existence of volatility fluctuations that are long-ranged correlated in 
time. Empirical studies show that the volatility correlation function decays very slowly in time, and is well fitted by 
a power-law of the time lag with a small exponent la la El Bl - An exponential fit with a single time scale - suggested 
by several simple stochastic volatility (SV) models |9 UlCl llU Il2l IT^ Il4 | - fails to describe the data: either the time 
scale is chosen to be small (a few days) and the long time tail is completely missed, or the time scale is chosen to be 
large (several hundreds of days), and the rapid initial decay is not properly accounted for. At least two time scales, 
and perhaps even three are needed to reproduce the volatility dynamics. Recent multifractal models 0, ITrj Il7| 
actually postulate the existence of an infinity of time scales to build a self-similar volatility model with a power-law 
autocorrelation function. 

Another, less well documented, effect is the negative return-volatility correlation: price drops are followed by an 
increased level of volatility. This effect was first noted by Black [18| and is usually called, somewhat inappropriatly, 
the 'leverage' effect, or asymmetric volatility correlation The presence of these correlations leads to skew in the 
distribution of returns and, correspondingly, to skew in the option smiles. The effect is much stronger for stock indices 
than for individual stocks, but appears to decay on rather short time scales (10-20 days ^3H3), much shorter than 
the long time scale noted above for the volatility correlation. 

Among the several SV models presented [Hi, both the Heston [ISIISO and thc Ornstcin-Uhlenbeck J§. 111 ItH] SV 
models are susceptible to contain an asymmetry parameter that generates some leverage correlations |21| . However, 
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as mentioned above, these models contain a single time scale which governs both the volatility correlations and the 
return- volatility correlations, and cannot account for the different temporal behaviour of the two. Multifractal models, 
on the other hand, have only been very recently extended to consider the leverage effect [zl,!^. These models are 
however quite complex and analytical results for the unconditional and conditional return distribution over different 
time scales are only partially available 22, '24| . 

Our aim in this paper is to extend the exactly soluble Ornstein-Uhlenbeck SV model introduced in 0, IT^, H"^ , 
such as to reproduce the multi-time scale dynamics reported above. In the present paper we allow for two different 
time scales, but the model can be generalised to account for as many time scales as needed, still retaining its exactly 
soluble nature. The basic idea is to make the mean reverting level of the volatility itself random and time dependent. 
The paper is divided into six sections. In section [n] we present our three-dimensional diffusion model. Sections lllll 
and IIVI respectively study the leverage and volatility correlations. In Section we apply the model to actual data, 
thus providing a way of estimating the parameters of the model. Conclusions are drawn in section IVll and some more 
technical details are relegated to appendices. 



II. A THREE-DIMENSIONAL DIFFUSION MODEL 



The dynamics of stock price S(t) is usually modelled as a multiplicative diffusion process. It is commonly assumed 
that stock price S(t) obeys the diffusion process: 

^.= f ,dt + adW 1 (t), (1) 

where /i is the drift, a is the so-called volatility and W%(t) is a Wiener process, i.e., dW\(t) — £i(t)dt, where £i(t) is 
Gaussian white noise with zero-mean and unit volatility: 

(£i(t))=0 and (6(*)6(0) = *(*-0- 
For the sake of simplicity we will work with the zero-mean return process defined by 

dX(t) e, |H - fxdt. (2) 

Market dynamics is then described by a simpler stochastic differential equation (SDE) 

dX(t) =adW x {t). (3) 

As is well-known, this model does not satisfactorily explain many properties of the market. A possible improvement 
is to assume that, besides the return following the diffusion process the volatility itself is a random variable, i.e., 
a = a(t) is a stochastic process. This is done in the stochastic volatility models which appeared in the literature 
during late eighties followin g si milar discrete time 'ARCH' models. There are several ways to choose the random 
process describing the volatility [l(J, [HI ■ The simplest one consists in assuming that the volatility is also a diffusion 
process following the Ornstein-Uhlenbeck (OU) process 0: 

da = -a(a-m)dt + kdW 2 (t), (4) 

where a > 0, k > 0, and m are nonrandom quantities, and dW2(t) is another Wiener process. The quantity 1/a is 
the characteristic time scale of the process and measures the strength of the reversion of <r(t) to the constant m which 
is the 'normal' level of volatility (see below). The quantity k is the amplitude of the volatility fluctuations and it is 
sometimes called 'vol of vol'. As is well known, the stationary distribution of the OU process is Gaussian and this 
implies that the probability density does not vanish for er — > 0. Since only |er| is meaningful, one has to 'fold' the 
negative part of the Gaussian onto the positive side with the result that the probability of | cr | has a small hump for 
| a | small. This and the Gaussian shape of the distribution of the volatility are not very realistic because the empirical 
distribution of the volatility is rather found to be close to a log-normal, or to an inverse Gamma distribution |25L l2q| . 
This is certainly a limitation of the present model. An interesting extension would be to consider the case where 
In cr, instead of a, follows an OU process. Another limitation of the present model is the assumption that dW\(t) is 
Gaussian. Although volatility fluctuations indeed generate fat tails, the empirical frequency of very large returns (like 
crashes) cannot be reproduced within the class of Gaussian stochastic volatility models, and some extra 'jumps' must 
be allowed (see the empirical results and comments in Section [VJ. 
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Two of us have recently studied the model represented by Eqs. JU-JU and showed 01 that it reproduces the 
leverage correlation only when the variations of volatility are anticorrelated with the variations of return. In other 
words, when both Wiener processes W\{t) and appearing in Eqs. @J and J3J are anticorrelated: 

(Hi(t)Ut'))=p5(t-ti), (5) 

where p < 0. As mentioned above the model given by Eqs. Q and J3J together with Eq. JSJ reproduces very well 
the observed leverage correlation but is unfortunately unable to reproduce simultaneously the long range volatility 
autocorrelation because leverage and volatility autocorrelation have both approximately the same characteristic time 
given by I /a. 

In the volatility process given by Eq. J3J, m is a nonrandom constant and indicates the level to which the average 
volatility eventually converges as time increases [27^ . In order to visualize this, suppose that the process begun in 
the infinite past (to — > —00) which implies that at time t the volatility has reached a stationary state, independent of 
initial conditions. In this situation one can easily see from the integration of Eq. J3J that 

a(t) = m + k [ e- a(t - t ' ) dW 2 {t'). (6) 



From this we clearly see why m can be called the normal level of volatility since (i) it is the stationary mean value of 
the volatility, (c(t)) = m, and (ii) m is the value to which a(t) would return without any random innovation (k = 0). 

However, there is no reason to believe that this 'normal' level of volatility is not itself slowly time evolving, with 
long periods where the average volatility is high (like stocks in the 30's), and other periods where it is smaller (like 
in the 60's). We will now relax this assumption and explore its consequences. We thus assume that m = m(t) is a 
random process also obeying an OU stochastic differential equation: 

dm = — a (m — m )dt + k^dW^{t), (7) 

where ao > and k$ > 0, and mo is now the 'true' long time average of the volatilty (but see the discussion 
in the conclusion). We therefore propose that the process describing the whole market dynamics is the following 
three-dimensional random process: 

dX[t) = a(t)dWi{t), (8) 
da(t) = -a[a(t) - m{t)]dt + kdW 2 (t), (9) 
dm(t) = -a [m(t) - m ]dt + k dW 3 (t), (10) 

where dWi(t) — £i(t)dt (i — 1,2,3) are Wiener processes, i.e., £,*(t) are zero-mean Gaussian white noises with cross- 
correlation given by 

mMt'))=piAt-t% (11) 

where pij = pji and pa — 1. We also suppose that Wz{t) is independent of W\(t) and W2(t). Hence, the correlation 
matrix is 



(/Hi) = \ P 1 (12) 




where — 1 < p < 1. As is common in finance, Eqs. J5|l- i|10|) are interpreted in the ltd sense and for the rest of the paper 
we will follow the Ito convention |2l|. A key consequence of the Ito convention is that <r(t) and m(t) are independent 
of any of the same time driving noises dWi(t) (i = 1, 2, 3). 

Note that the choice of the correlation matrix given by Eq. I|12l) assumes that the variations of return and instan- 
taneous volatility are directly correlated while the variation of the volatility normal level is unaffected by recent price 
changes. Although reasonable from a financial point of view, there exists a different limit in which the normal level of 
volatility is directly correlated with the return. This case is treated in Appendix El where we focus on an alternative 
correlation matrix for which pi 2 — and P13 = p instead of p\ 2 = p and P13 = 0. We will see there that in this limit 
both the leverage and volatility correlations are governed, on long time scales, by the 'slow mode' of the model. 

In what follows we will assume that at time t the entire process has reached the stationary state. From Eq. JJJ we 
see that the expression for m(t) in the stationary state is 

m(t) = m + k f e-^-^dW^t'). (13) 
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Now the stationary expression for m(t) is given by the solution of Eq. © rather than by Eq. (|BJ). We thus have 

a{t)= f \k&{t') + amtf)]e- a ( t - t ">dt, 

J — oo 

which, after using Eq. I|13|> . reads 



a(t)=m + k e- a{t ' t ' ) dW 2 (t') + 



kg 
1 - A 



e -a (t-t') _ e -a(t-t') 



dW 3 (t'), 



where 



A = — . 

a 



(14) 



(15) 



We will see later on that empirical evidence implies that < A < 1, as suggested by the initial motivation of the 
model. One of the interest of this linear model is that one can also compute explicitely, using its Markovian nature 
the time evolution of the volatility conditioned to a given present value of the instantaneous volatility (J (to) and of 
the 'normal' level of the volatility m(t ). This in turn opens the way to obtain the conditional distribution of returns 
over an arbitrary horizon T, a quantity needed for pricing options (see [9l Il2l IT^V 
Observe that since m(t) is a Gaussian process its mean and variance read 



and the autocorrelation is 



(m(t)) = m , Var[m(i)] = in^v 2 , 
(m(t)m(t + t)) = ml (l + vle~ a ° T ) (r > 0), 



where r > and v§ is defined by 



2rriQaQ 



Likewise a(t) is also Gaussian and has the following averages 



(a(t)) = m , Var[(r(t)] = m 2 [v 4 + v& 



and autocorrelation 



where t > 0, and 



(a(t)a(t + r)) = m 2 



y 2 _ H 



1 - A 2 



1 + A 



1 - A 2 



(16) 



(17) 



(18) 



(19) 



(20) 



2TOqQ! 



(21) 



III. THE LEVERAGE EFFECT 

As recalled in the introduction, it has been known for long that volatility and returns are negatively correlated 
|l8j |. This correlation has however not been studied from a quantitative manner until quite recently |2(1 |. where the 
following 'leverage' correlation function was studied: 



C(r)^^(dX(t + T) 2 dX(t)) 



(22) 



where dX(t) is the zero- mean return and 



2\ 2 



Z = (dX(t) 2 } 



(23) 
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is a convenient normalization coefficient. A large amount of daily returns for both market indices and share prices 
were studied in |20|. with the finding that, 

(A,T > 0). Hence, there is an exponentially decaying negative correlation between future volatility and past returns 
changes. No correlation is found between past volatility and future price changes (except perhaps for very short time 
lags |r|). In this way, a clear causality of the leverage effect is established, whereas there are conflicting claims on 
this issue in the literature For individual stocks, the coefficient A is found to be close to 2, a value expected on 
the basis of simple 'retarded' effect, where the price changes are locally additive, rather than multiplicative random 
variables |2jj. For indices, however, A is much stronger (A w 20), and the retarded interpretation fails. The coefficient 
r, on the other hand, defines the inverse relaxation time of the leverage correlation, and is found to be on the order 
of ten days for stock indices. 

Let us study the leverage correlation within the framework of the three-dimensional diffusion model presented in 
Section [H] First, from Eq. ©, we calculate the correlation 

(dX(t)dX(t + t) 2 ) = {a{t)dWi{t)a{t + rfdW^t + r) 2 ). (25) 

Following the Ito convention, if r > then dWi(t + r) is uncorrelated with the rest of stochastic variables. Thus, 
taking into account that 

(dWi(t + r) 2 ) = dt, 

we have 



(dX(t)dX(t + r) 2 ) = (cr(t)dWi (t)a{t + r) 2 )dt (r > 0). 

Otherwise, if r < 0, then dW\{t) is uncorrelated with the remaining variables. Hence, taking into account that 
(dWi(t)) = 0, we get 

(dX(t)dX(t + r) 2 ) = (t<0). 

Therefore, 

(dX(t)dX(t + t) 2 ) = (a(t)a(t + T) 2 dW!(t))R(T)dt, (26) 

where 

h m={J; < 27 > 

is the Heaviside step function. Now the Novikov theorem allows us to write (see Appendix |A"|) 

(dX(t)dX(t + t) 2 ) = 2pke- aT R(T)(a(t)a(t + r))dt 2 , (28) 

where (a(t)a(t + r)) is given by Eq. 1(217)1. 

On the other hand, as a consequence of the Ito convention, processes a(t) and dW\(t) are independent. Therefore 
from Eqs. JSJ) and ill' Ml we see that 

Z = (dX(t) 2 ) 2 = (aitffidWxit) 2 ) 2 = (a(t) 2 ) 2 dt 2 . 

Taking into account Eq. (flfjf) . Z reads 



Z = mUl + S + ^L_) dt 2 , (29) 

Finally, the substitution of Eqs. (28) and ^ into Eq. JUJ yields 

C{T)=-H(T)A( T )e- aT , (30) 
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where 



with: 



and 



A{r) e, 2pvV ^ [1 + oe— + fee—] , (31) 

m (1 + a + b) 



01 = (33) 

Note that Eq. 1|3()[) is not very different from Eq. (|24(l above: the dominant time decay is fixed by the short time scale 
of the volatility process 1/a. As intuitively expected, the correlation coefficient p between dW\ (returns) and dWi 
(volatility changes) needs to be negative if we want to describe the observed anticorrelation. 

IV. VOLATILITY CORRELATION 

The correlation for the volatility a(t) has been already given by Eq. i|2(J|) . Nevertheless, in order to study the 
different time regimes appearing in the correlation, it is convenient to deal with a slightly different form of correlation. 
We will thus work with the function 

CorxK*), a(t + r)]= «'(*M' + ')» 



VVar[<7(t)yVar[<r(t + r)] 
where ((a(t)a(t + r))) is the second cumulant of the volatility process defined by 

({a(t)a(t + r)» = (a(t)a(t + r)) (a(t)){a(t + r)>. 

From Eqs ifT^ll and (PU|l we have 

Corr [a(t),a(t + r)] = M[ae~ aT + 6e~ Q ° T ], (34) 

where 



M = ~r^, (35) 



and Vq is given by Eq. 

If we assume that a ^ <* the volatility correlation (|34|l has two different regimes according to r being small or 
large. 

• Suppose that r is small and such that the following two conditions meet simultaneously 

aoT <C 1 and ar ~ 1. 

In this case, the characteristic time scale of volatility process given by 1/a (see Eq. ©) will dominate and the 
correlation reads 

Corr [o-(t), a{t + t)] - M(b + ae- aT ), (a r < 1). (36) 

• Suppose now that r is large and such that 

ar 3> 1 and o.qt ~ 1. 

In this case the characteristic time of the normal level of volatility m(t) (see Eq. (|10f) ') is the dominant one and 
we see from Eq. (|34ll that 

Corr [a{t], a{t + r)] - Mbe~ a ° T , (ar > 1). (37) 
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Note that both approximations are consistent since we have assumed that «o < a. This implies that A <C 1 (cf. 
Eq (jlSp) is very small which, in turn, simplifies the form of the parameters M, a, and b. Indeed, now we can write 

M= i^ + u^)- 1 +0(A), (38) 

and 

a = v 2 + 0(A), b = t% + 0(\). (39) 

Equations l|36(l and l|37|) show respectively the correlation behavior for small and large time r. The long range 
correlation for the volatility is governed by ao, that is, as expected, by the relaxation process of the volatility to its 
long time average uiq. This behavior is different from the leverage correlation which has a faster time decay a -1 that 
sets the 'fast' reversion of the volatility. Notice that this is consistent with empirical observations |^ll5ll27| (see also 
Fig. n and discussion in the next section) . 

At this stage, we should remind that the actual evaluation of the volatility a(t) is very difficult, since the volatility 
itself is not observed. In practice one can measure a noisy measure of the instantaneous volatility from: 



\<j(t)\^^J[X(t + At)-X(t)} 2 /At, (40) 

where X(t) is the (undrifted) log-price defined above. Therefore, only \<r(t)\ can be extracted from data, and this is in 
principle insufficient to estimate (|34|l . Another possible expression for the empirical autocorrelation of the volatility 
that is more directly observable is: 

Corr [dX{t) 2 ,dX{t + r) 2 ] = y((dX(t) 2 dX(t + t) 2 )), (41) 

where Y 2 = Var[dA(i) 2 ]Var[dA(i + r) 2 ] and 

((dX(t) 2 dX{t + T) 2 )) = (dX(t) 2 dX(t + r) 2 ) - (dX(t) 2 )(dX(t + T) 2 }, 

is the second cumulant associated with dX(t) 2 . This correlation function is particularly interesting because it is 
directly related to the kurtosis of the terminal distribution of price changes, a quantity which is in turn related to the 
curvature of volatility smiles [2^| . 

In Appendix [5] wc show that within the present Gaussian model, the empirical correlation coefficient given by 
Eq. I|41|) can be written in terms of the volatility autocorrelation as 



and from Eqs. I|19|l and (|20[1 we explicitly write 



a(2 + ae- aT )e~ aT + b(2 + be- a ° T )e~ a ° T + 2abe- {a+a ° )T 



(43) 



Corr [dX(t) 2 ,dX(t + r) 2 ] = N 

where a and b are given by Eq. (|32[) and 

N = [1 + 8(v 2 + D 2 ) + 4(*/ 2 + D 2 ) 2 ] ^ . (44) 

Observe that correlation (|43|) consists of a linear combination of decaying exponentials with arguments (XqT, ar, 2aoT, 
2ar, and (ao + a)r. We study some asymptotic limits for the correlation (f^ as done above with a{t). We again 
consider the case when ao <C a {i.e., A <C 1) and study the short and long r behavior. For small t, when aoT is small 
but ar is finite and not too small, we have 

Corr [dX(t) 2 ,dX(t + r) 2 ] ^ N [(2 + b)b + a (2 + 2b + ae~ aT ) e~ QT ] , (a r < 1). (45) 

In this regime, we only have exponentials with a as it happened with the sigma autocorrelation (cf. Eq. I|36l) ). For 
large t, the correlation is again dominated by the exponential with characteristic time scale ao and reads 

Corr [dX(t) 2 ,dX(t + r) 2 ] w Nb (2 + b e - a " T ) e~ Q ° r (ar > 1). (46) 

Observe that pairs of Eqs. (|36|l - (|37|l and Eqs. I|45|l - (|46|l contain the behavior that is empirically observed in financial 
markets. Although Eq. (|43|l has a much more complicated form, it still contains the two desired regimes for large and 
small r as we have discussed above. 




time (in days) 

FIG. 1: The autocorrelation of AX(t) 2 for the Dow- Jones daily index. We also plot their theoretical correlations given by 
Eq. 1431 with the parameters estimated from the fit with a = 0.1 days -1 , ao — 1.3 x 10~ 3 days -1 , a = 0.14 and b — 0.04. The 
inset is a zoom of the correlation focussing on the range from 1 to 80 days. 



V. ESTIMATION OF PARAMETERS 



Any acceptable market model must be easily calibrated to reproduce the dynamics of the market. In other words, 
the model should be supplemented by a systematic methodology for estimating its parameters. In our model we have 
six parameters to estimate. These parameters are p, k and fco (or v 2 and Vq), a and ao, and too (see Eqs. i[%|l-i[10|l). 

Let us discuss the way to estimate these parameters. We will proceed in a very similar way to the case given 
in Ref. |13| and perform the empirical estimation of the Dow-Jones index daily returns by approximating dX by 
AX = X(t + At) — X(t), i.e., 

dX(t) w X(t + At) -X(t), 

where At = 1 day. From Eqs. I|19[) and l|B9|l and taking into account that a(t) is independent of dWi(t), and that 
dW\(t) is Gaussian, we have 

Var[AA] = {a{t) 2 )At = ml{l + v 2 + vl)At 1 (47) 
Var[AX 2 ] = [3(a(t) 4 ) - (a{t) 2 )] At 2 =2m% [4(1 + ^ 2 + *> 2 ) 2 - 3] At 2 , 

where v 2 = fc 2 /(2am§) and v 2 = k 2 /(2aorn^{l + A)). From these equations we get 

1 _ 4 1 Var[AX 2 ] 

{l + v 2 + vl) 2 3 6Var[AX] 2 ' ^ ' 

Hence, we can estimate the values of (v 2 + v 2 ) once we know the empirical values of the daily variances of A A" 
and AX 2 . In a second step, mo is estimated from the knowledge of [y 2 + i> 2 ) and the empirical value of Var[AA]. 
In Table [IJ we briefly report these operations and give the corresponding estimation of iy 2 + v 2 ) and Too for the 
Dow- Jones index time-series from 1900 until 2000. 
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FIG. 2: The leverage effect in the Dow- Jones daily index. We plot the leverage function C(t) for the Dow- Jones index 
from 1900 until 2000. We also plot a fit with the leverage function H30H with a — 0.1 days -1 , a = 0.14, b — 0.04, and 
m = 1.19 x 10~ 2 days~ 1/2 finding that p = -0.48 (see Table HD- 

TABLE I: Estimate from return variances. We estimate parameters of the model from Dow- Jones historical daily returns from 
1900 to 2000. We use the identity 14811 and take the variance 1471 to derive the estimated quantities (y 2 + u 2 ) and mo. 



Estimators Dow-Joncs daily return Theoretical values 



Var[AX(t)] 
Var [AX(tf] 



1.68 x 10 -4 
10.5 x 10" 8 



2mt[^u 2 + u 2 ) 2 



3] At 2 



Parameter estimates 
.2 1 ~.2 



mo 



0.18 
18.9% year" 172 



Since (v 2 +t%) is always positive and (1 + v 2 + v 2 ) 2 > 1, we see from Eq. (|48f) that 



and the kurtosis defined by 



Var[AX 2 ] 
2 " Var[AX] 2 < 8 ' 

_ Var [AX 2 ] 
K ~ Var[AX] 2 ~ 2 
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TABLE II: Parameter estimation calculated from the volatility autocorrelation (1431) and the leverage correlation (1301) . We 
estimate a, z>q, and v 2 from the empirical volatility autocorrelation of Dow- Jones stock index plotted in Fig. Q We estimate 
the parameters p and check a from the fit of the Dow- Jones stock index data leverage plotted in Fig. [5] 



Estimators Dow-Jones data estimation 

Volatility autocorrelation 
a 0.1 days -1 

q 1.3 x 10" 3 days" 1 

a~u 2 0.14 
6 ~ i/g 0.04 

\ = a /a 1.3 x 10 -2 

k = m V2av 2 2.0 x 10~ 3 days -1 

k = m y / 2Qoi'o 1-2 x 10 -4 days -1 

Leverage 

£(0+) -11.9 
p -0.48 



has the bounds 

< k < 6. (49) 

As expected, volatility fluctuations induce some excess kurtosis, but this excess kurtosis cannot take arbitrarily high 
values. For the Dow-Jones index data and after taking away 'crash days' of amplitude greater than 12%, we get a 
kurtosis of k ~ 1.7 which is consistent with the requirement (|49|l . The removal of these outliers only represents 11 
data points, a very small fraction of the total amount (more than 27000 sample data points). However, as many 
authors have reported H^,|33,|^|, the total kurtosis is found to be larger, perhaps even infinite (in particular for 
higher - intraday - frequencies) . This is a limitation of the present model for describing higher frequency movements 
where one should eventually take into account a non Gaussian character for dW\ (see the discussion in |26) and 
Refs. IHIll) 

We have four more parameters to estimate, a, ao, v 2 and p, and we will use both the leverage and the volatility 
correlations to obtain them. As we see in Figs. ^ and El the leverage correlation appears to be much more noisy than 
the volatility autocorrelation. For this reason we start with volatility autocorrelation for quantifying a, cvo, v 2 . We 
subsequently analyze the leverage correlation for testing the estimated a and for finally obtaining the parameter p. 

As we have shown in the Section llVI there are different ways to interpret the volatility correlation from empirical 
data. In the previous Section our choice was to take AX(t) 2 as the square of the daily volatility. We have computed 
the empirical volatility correlation and plotted it in Fig. ^ with the corresponding fit with our theoretical correlation 
in Eq. g3J): 



Corr [dX(t) 2 , dX{t + t) 2 ] = N 



6(2 + be- aoT )e- aaT + 2abe' [a+ao)T 



Ai>§/(1 - A), b = t%/(l - A), N = [1 + 8{v 2 + i> 2 ) + 4(^ 2 + £ 2 ) 2 ] . The fit plotted in Fig. Q] is done 
with the constraint that a + b = (y 2 + i> 2 ) = 0.18. 

The proposed fit is obtained by taking l/an — 3 years and l/a ~ 10 days. Recall that A = ao/a is been assumed 
in previous Sections to be very small and this is now confirmed by empirical observations since A = 1.3 x 10~ 2 . In 
such a case, we can consider that a ~ v 2 and b ~ fg. The results obtained are given in Table HP jointly with k and ko 
that are given by Eqs. (JTSJ and J2U- 

The leverage effect provides a way to estimate the correlation coefficient p. It also serves us as a test of consistency 
in the estimation of the characteristic time l/a. In a previous work [l3|, it is been suggested that l/a ~ 20 days. 
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FIG. 3: Dow- Jones index and simulation paths over 10000 days which corresponds to approximately 40 years. Top figure 
shows the Dow- Jones index daily return changes over 10000 trading days and compare it with several simulations. From top 
to bottom we show random path of our three-dimensional stochstic volatility model (ISt- dlOIL the path of the two-dimensional 
approach and the Wiener process assuming constant a. Parameters are given in Tables HI and HT1 

However, the fit of the leverage function (|3U|) — (|31 |> to the Dow- Jones daily data (Fig. [2J is not very binding. The 
quality of this non linear fit remains basically unchanged for modifications in l/a inside the range 7 — 20 days. For 
this reason we will keep our estimated a from volatility autocorrelation and accept that the results from the leverage 
correlation are consistent with the ones derived from the volatility correlation. 

Let us now focus on the estimation of p. Recall that our theoretical leverage is given by Eqs. (|3U|1 - (|3 1|) with a and 
b given by Eq. (|32|) . When r — > + , we have 



£(0+) 




(50) 



(1 + v 2 + v^)m ' 



which allows us to estimate p once we know a = 0.1 days 1 and (y 1 + v$) = 0.18. All these operations are summarized 
in Table ITT1 The estimated magnitudes C(0 + ) and a obtained from the Dow-Jones index are of the same order as 
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FIG. 4: Dow- Jones index and simulation paths over 1000 days (approximately 4 years). Top figure shows the Dow- Jones index 
daily return changes over 1000 trading days and compare it with several simulations. From top to bottom we show random 
path of our three-dimensional stochstic volatility model (ISIl- llOII . the path of the two-dimensional approach and the 

Wiener process (J^J assuming constant a. Parameters are given in Tables HI and HTl 



those given in [2(| from an average over several stock indices. The value of a is also consistent with that found 
by Dragulescu and Yakovenko within the Heston model 0] (for which the volatility correlation function is a single 
exponential) . 

Finally, in Figs.|21and2|we simulate the resulting process with the parameters estimated above. Once we know the 
six estimated parameters, we can simulate the random dynamics for AX(t) and compare it with the empirical Dow- 
Jones time series. We also simulate the ordinary geometric Brownian motion where the volatility is constant together 
with the two-dimensional SV model studied in [13j in which the normal level of volatility is constant. In Figs.[3]and0] 
we see that the ordinary geometric Brownian model cannot describe both the largest and the smallest fluctuations of 
daily returns. Furthermore, the two-dimensional SV model is not able to describe the clustering that appears when 
we take a longer number of days (see Fig.[3J). Hence, we can conclude that SV models generate trajectories similar to 
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that of the Dow- Jones but the three-dimensional model presented herein improves this description when we observe 
a wider time window. On the other hand, one should keep in mind that some aspects of the empirical data, such as 
the nearly log-normal distribution of volatility, or the extreme tails of the distribution of returns, are not reproduced 
by the present version of the model. 

VI. CONCLUSIONS 

In this paper, we have suggested a simple way to build a stochastic volatility model that accounts for the presence of 
multiple time scales in financial markets. The volatility correlation function contains at least a 'short' time scale (tens 
of days) and a very long time scale that does not seem to be present in the asymmetric return-volatility, which only 
exhibits the 'short' time scale. The idea is to let the 'normal' level of volatility be itself time dependent and evolve 
over long time scales. The simplest setting is that of the linear Ornstein-Uhlenbeck model, where the calculations are 

simple mm mi, 

and which can easily be calibrated to market data. The overall description of the volatility and the 
leverage correlation using this three dimensional model is quite satisfactory, although several deficiencies remain: (a) 
the volatility process appears to need at least a third, very short time scale (on the order of a day), to be satisfactorily 
accounted for, (b) the empirical distribution of volatility is poorly reproduced by the present Gaussian model, and 
(c) the assumption that dW\ is Gaussian does not allow one to account for extreme events. There are several natural 
extensions of the present model that could deal with these discrepancies: (i) one could generalize the Heston model as 
is done here, such as to avoid the unrealistic 'hump' of the volatility distribution close to zero; (ii) one could consider 
a model where the log- volatility, rather than the volatility, obeys an Ornstein-Uhlenbeck process, and (iii) one could 
build a 'Matrioshka doll' Ornstein-Uhlenbeck process where the construction given here is iterated as many times as 
needed to reproduce all the empirical time scales. In the limit of an infinite number of time scales, it might be possible 
to construct in an alternative way the multifractal random walk of Bacry, Muzy and Delour. Finally, it would be 
useful to extend the computation of the present paper to compute the conditional distribution of price changes over 
different time scales, such as to obtain option prices. We hope to report the results of this computation in the near 
future. 
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APPENDIX A: THE NOVIKOV THEOREM 

The Novikov theorem provides a way of evaluating averages involving the Wiener process and functionals of it. 
Suppose we have a well-behaved function of time f(t; [£]) also depending functionally on an rt-dimensional Gaussian 
white noise £(t) = • ■ ■ ,£„(£)). The Novikov theorem (see Ref. |3J] for further details) states that the nonlinear 

average {fit; can be written as 

</(*; [£])&(*')> = (^r0y) > Cj = 1,2,3,. ■ ■), (Ai) 

where 8f(t; [$])/5^j(t') stands for the functional derivative of / with respect to £ j^]. We can write an alternative, 
and sometimes more convenient form of the theorem if we specify the correlation between the n components of 

&(ti)t j (t2)) = p ij 6(t 1 -t 2 ). 

One can easily see that in this case we can write 

&(*)=£(*) (A2) 

where £(t) is a Gaussian white noise independent of any other component of £(t), i.e., (^(t)^j(t')) = for any j ^ i 
and t' . On the other hand the chain rule applied to the functional derivative allows us to write [35| 
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Then, taking into account Eq. (|A2|) . the Novikov theorem can be set in the following alternative form: 



(/(*;[£])&•(*')> = £ 



/<*/(*; [€]) 



(A3) 



We are now able to prove Eq. (|28[1 . In this case we have to consider the average (a(t)a(t + t) 2 £i (£)). Consequently 
[£]) = &(t)a(t + r) 2 and the Novikov theorem, in the form given by Eq. IjAlfl . tell us that 



On the other hand we know from Eq. (|14() that 



S [a(t)a(t + t) 2 



a(t) = mo + fc / ^(t , )e- a ( t -*'W + 



fen 



1-A 



&(f) 



e -ao(t-*') _ g-a(t-t') 



dt'. 



(A4) 



(A5) 



where the integrals are represented in the Ito sense. Observe that <j(t) depends functionally on £a(t) and ^(i) but not 
on and, as it is shown by Eq. Q12[l. and ^(i) are the only two correlated process with correlation coefficient 
p. Thus, the functional derivative reads 

5 [a(t)a(t + rf] 5 [cr(t)cr(t + r) 2 ] 
— = p- 



Sti(t) r %(t) 

where we have used Eq. (|A2(I to see that d^2(t)/d^(t)i = p. Therefore, 



6 [a(t)a(t + r) 2 



= P 



(A6) 



The functional derivatives of a(t) and cr(f + r) with respect to ^(i) can be obtained through Eq. (|A5fl after taking 
into account that ^t(t) is independent of ^(t) and that d^2(t')/S^2(t) — S(t — t') j^. Thus 



fcH(O) = 0, 



(A7) 



since H(0) — (cf. Eq. 1271 1. Note that this is consistent with the Ito interpretation which assumes that ^s(t) lies 
outside of the integration interval (— oo, t) of the first integral on the rhs of Eq. I|A5() . Analogously 



6a(t- 



fce- QT H(r). 



(A8) 



Finally, from Eq l|A6[) . we have 



(a(t)a(t + r)%(t)) = 2pke- aT (a(t)a(t + r))H(r), 



which proves Eq. I)28[l. 



APPENDIX B: DERIVATION OF EQUATION (02) 

We will demonstrate the equality 

(dX(t) 2 dX(t + t) 2 ) - (dX(t) 2 ) (dX(t + t) 2 ) _ (a(t)a(t + r)) 2 - (trft)) 4 
^Var[dX(<) 2 ]^Var[rfA(i + t) 2 ] ~ 4(a(i) 2 ) 2 - 3(ct(0) 4 



(Bl) 



which is equivalent to prove Eq. I|42|l from Eq. (|41|l . To this end we will deal with each term on the left hand side of 
Eq. HB1|) separately. 

We first derive (dX(t) 2 dX(t + r) 2 ). From Eq. © we have 

(dX(t) 2 dX(t + T) 2 ) = (aitfait + TfdW^t^itydt 2 , (B2) 
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where we have assumed that dWi(t+r) is independent of the rest of stochastic variables and also that (dWi(t + T) 2 ) — 
dt. Since dW\{t) — £,\{t)dt we can write the average on the rhs of Eq. (IB2[> as (a(t) 2 a(t + r) 2 : dWi(t)£i(t)) . Now the 
Novikov theorem, Eqs. IjAlj) and (|A3|) . yields 



(a(t) 2 a(t + T) 2 dW 1 (t)£ 1 {t))=2p 



5a{t) a(t)a(t + r) 2 dWi{t)\ + / a(t) 2 a(t + r ) ^±fW (t) 



6&{t) v v ■ v 1V 7 \ w %(t) 

+ (cr(t) 2 cr(t + r) 



2 , va*[<Wi(*)] 



«i(t) 



which, after using Eqs. (|A7(I - (|A8|I . reads 

(a(t) 2 a(t + r) 2 dW 1 m 1 (t))=2pke- aT R(r)(a(t) 2 a(t + r)dWx (t)) + (a(t) 2 a(t + r) 2 ^ffl ) ■ 

In order to obtain an expression for the functional derivative 5[dWi(t)]/5£i (i), we write dW\ (t) in a somewhat intricate 
form 

t+dt 



^(t')dt' 



Thus 
Hence 

(ct(£) 2 ct(£ + r) 2 dlFi(t)Ci(t)) = 2 / 9fc e - QT H(r)(CT(<)CT(t + r) 2 dWi(*)) + (<K*)M* + ^) 2 ). (B3) 
We apply again the Novikov theorem to the average (a(t)a(t + r) 2 dW\(t)) with the result 

(a(t)a(t + r) 2 dWi(t)) = p I -/j-r[a(t) 2 a(t + r)]\ dt, 



56 (*)' 

which, after using Eqs. l|A"7 )) -l|A"8 |l . yields 

(a(t) 2 a(t + T)dWi(t)) = pke- aT R(T)(a{t) 2 )dt. 

Substituting this into Eq. 1|B3(1 we have 

(a(t) 2 a(t + r) 2 dWi(t)&(i)) - 2p 2 k 2 e- 2aT R(r)(a(t) 2 )dt + {a(t) 2 a{t + r) 2 ). 

Notice that due to the differential dt, the first term on the rhs of this equation is negligible in front of the second 
term. Thus 

(a(t) 2 a{t + T) 2 dWi(t)^(t)) = (a(t) 2 a{t + r) 2 ) + 0(dt), (B4) 
and the substitution of Eq. (|B4|I into Eq. (|B2() yields 

(dX(t) 2 dX(t + t) 2 ) = [{a{t) 2 a{t + t) 2 ) + O(dt)] dt 2 . (B5) 
On the other hand, since dW\(t) is independent of a(t) we can write 

(dX{t) 2 ) (dX(t + t) 2 ) = (<r(t) 2 ) (a(t + T) 2 )dt 2 
and, taking into account that process is in the stationary regime, we have 

(dX(t) 2 )(dX(t + t) 2 ) = (a(t) 2 ) 2 dt 2 . (B6) 

Finally, we derive the variance: 

Var[dX(t) 2 ] = (a{t) 4 )(dW(t) 4 ) - (a(t) 2 ) 2 (dW(t) 2 ) 2 , 
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where we take into account that a(t) and dW\(t) are independent. We can go further since (dW(t)' 1 ) — 3(dW(t) 2 ) 2 = 
3dt 2 . Hence, 

Var[dA(<) 2 ] = 3(a{t) 4 )dt 2 - (a(t) 2 ) 2 dt 2 , (B7) 

and, due to the fact we are dealing with process in the stationary regime, the same result applies to the variance 
Var[dX(t + r) 2 ]. 

Therefore, from Eqs. I|B5(1 . (|B6|I and (|B7() . we have that correlation function given by Eq. CE) is 

Corr [dX (t) 2 dX(t + rf] = ^ff'^ + (*)- ( B8 ) 

This expression contains two terms which must be also derived: (a(t) 2 a(t + r) 2 ) and (<j(t) i ). Again, these terms can 
be obtained with the help of Novikov theorem. Lengthy and tedious calculations allow us write them in terms of 
other expressions presented above (c/. Eqs. H19|) and (J20J|). The first one is 

{a(t) 2 a(t + r) 2 ) = 2 [{a(t)a(t + r)) 2 (a(£)) 4 ] + (a(t) 2 ) 2 , 

while the second one is obtained from the first one by assuming r = 0. That is: 

{o{tf)=Z{o{t) 2 ) 2 -2{a{t))\ (B9) 

If we insert these averages to the correlation (|B8|) . we prove the equality given by Eq. (|B1|) and finally obtain Eq. lH2"|) . 



APPENDIX C: AN ALTERNATIVE NOISE CORRELATION 



We recall that our model is given by a three-dimensional random process © (I1UII 

dX{t) = a(t)dWi(t), 

da(t) = -a[a(t) ~ m{t)]dt + kdW 2 (t), 

dm(t) = — ao[m(t) — mo]dt + kodW^it), 

where dWi(t) = £,i(t)dt (i = 1,2,3) are Wiener processes, i.e., ^(t) are zero-mean Gaussian white noises with cross- 
correlation given by (^i(t)^j(t')} = pijS(t — t'). Note that pij — pji and pa — 1. We have supposed that Wz{t) is 
independent of W\ (t) and Wz (t) . And we also have assumed that there exists a non zero correlation between W\ (t) 
and Wz(t) that was quantified by the parameter p. 

This appendix modifies this framework by changing the existing correlations betweeen the Wiener processes Wi 
(i = 1, 2, 3). We want to analyze the consequences on the model when we allow correlations between W\(t) and W^(t) 
instead of between Wi(t) and Wzit)- In this case, the correlation matrix is 




(Pij)= 10 (CI) 



where, again, — 1 < p < 1. 

One can easily see that this change in the correlation matrix does not affect most of the main statistical properties 
of the model although there is one important property that is modified: the leverage correlation. 
According to the definition 1221) , the leverage correlation is quantified by 

£(t) = ]-{dX{t + T) 2 dX{t)) (C2) 

where X(t) is the zero-mean return and Z = (dX{t) 2 ) 2 = mg(l + v 2 + vffidi 1 (cf. Eq. (|2~9"|> ). We need to calculate 
the correlation 

(dX(t + TfdX{t)) = H{T){(T{t)dWl{t)(j{t + TfdW X (t + T) 2 ). 

Its derivation is similar to the one given in Section fllll Since we follow the ltd convention, dW\{t + r) is independent 
of the rest of variables when r > 0. This fact and the averages (dW\(t)) — and (dWi(t) 2 ) = dt allow us to show 
that 



(dX(t + T) 2 dX{t)) = (a(t)a(t + t) 2 dW^H^dt, 




time (in days) 



FIG. 5: The leverage effect in the Dow- Jones daily index with the alternative leverage. We plot the leverage function £(r) and 
£(t) and compare with the Dow- Jones index from 1900 until 2000. The fit with the leverage function C(t) given by Eq. JSUJl 
takes parameters given by Table IH1 We also plot the two alternative leverage functions given by Eq. fl(J7ft . Alternative 1 shows 
the possible leverage with a — 0.1, a = 0.14, Qo = 0.04 and b = 0.04 while alternative 2 takes parameters qq = 0.1, b = 0.14, 
a = 0.04 and a = 0.04. 



where H(t) is the Heaviside step function given by Eq. (|27fl. The next step of the calculation needs to apply the 
Novikov theorem enunciated in Appendix [S] although now its result appears to be different from the one given by 
Eq. H28|) since the correlation is between W\{t) and W%(t), In the present case the Novikov theorem given by Eq. (|A3|) 
allows us to write 



{a(t)a{t + T) 2 dW 1 (t))=J2pu 

i=i 

On the other hand we know from Eq. (|14l) that 



S[a(t)a(t + T) 2 } 



dt. 



a(t) =m + k / ^(t')e- a ^dt' + 



kg 

1 - A 



6(f) 



e -ao(t-t') g-a(t-t') 



dt' 



(C3) 



(C4) 



where the integrals are represented in the ltd sense. Since volatility (|C4|) does not depend on £i and the only two non 
zero contributions correspond to p\\ = 1 and p\z = p (c.f. Eq. I|C1||1. 



But from Eq. Id 



(a(t)a(t + T ydWi(t)) = 2p { a(t) 



6a(t + t) fco 



1 - A 



t 5a{t + T) 



H(r), 



dt. 



(C5) 



(C6) 
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and (cr(t)) = toq- Hence the leverage correlation given by Eq. I|C2|) reads 



C(r) =H(r)A(r)(e" 



where 



A(r) 



m (l 



A)(l + ^ + £>2) 2 



3±_ 

1 - A 



1 - A 



(C7) 



(C8) 



We now compare the two leverage correlations C(t) and C(t). In Section ITTT1 we obtained C(t) given by Eqs. 
(EU) and see that the limit £(0+) goes to a certain negative value. In contrast, in the present alternative model, 
£(0 + ) = 0. In Fig- El w e see that leverage of the Dow- Jones index goes to —11.9 when t — > + . Hence, this alternative 
model is in clear contradiction with empirical observations. In addition, we also observe in Fig.0that £(r) has a longer 
range than the one observed in the Dow-Jones index. The reason for this being is that the exponential time decay of 
the leverage function l|C8|) is dominated by the same parameter as to the case of the volatility autocorrelation 143fl . 
Thus, in the present case, leverage and volatility correlations have the same long range persistence and this is also in 
clear contradiction with the empirical results telling us that leverage has a much shorter memory than the volatility 
(see and compare Figs. Q and [SJ. 



[1] M. Dacorogna, R. Gengay, U. Miiller, R. Olsen, and O. Pictet, An Introduction to High- Frequency Finance (Academic 
Press, London, 2001). 

[2] R. Cont, Empirical properties of asset returns: stylized facts and statistical issues, Quantitative Finance 1 223 (2001). 
[3] D. M. Guillaume, M. M. Dacorogna, R. D. Dave, U. A. Miiller, R. B. Olsen and O. V. Pictet, From the bird's eye to the 
microscope: A survey of new stylized facts of the intra- daily foreign exchange markets Finance and Stochastics 1 95 (1997). 
[4] R. Mantegna, H. E. Stanley, An Introduction to Econophysics (Cambridge University Press, Cambridge, 1999). 
[5] A. Lo, Long term memory in stock market prices, Econometrica 59 1279 (1991). 

[6] Z. Ding, C. W. J. Granger, and R. F. Engle, A long memory property of stock market returns and a new model, Journal 

of Empirical Finance, 1 83 (1993). 
[7] Y. Liu, P. Cizeau, M. Meyer, C.-K. Peng, H. E. Stanley, Correlations in Economic Time Series, Physica A 245 437 (1997). 
[8] J. F. Muzy, J. Delour, and E. Bacry, Modelling fluctuations of financial times series: from cascade process to stochastic 

volatility model, The European Physical Journal B 17 537 (2000); E. Bacry, J. Delour, J. F. Muzy, Multifractal random 

walk, Physical Review E 64 026103 (2001). 
[9] E. Stein and J. Stein, Stock Price Distributions with Stochastic Volatility: An Analytic Approach, Review of Financial 

Studies 4 727 (1991); reprinted in Volatility: New Estimation Techniques for Pricing Derivatives, Robert Jarrow Ed., Risk 

Publications, 1998, pp. 325-339. 
[10] S. Heston, A closed-form solution for options with stochastic volatiliy with applications to bond and currency options, 

Review of Financial Studies 6 327 (1993). 
[11] J. -P. Fouque, G. Papanicolaou, and K. R. Sircar, Derivatives in Financial Markets with Stochastic Volatility (Cambridge 

University Press, Cambridge, 2000). 
[12] R. Schobel and J. Zhu, Stochastic volatility with an Ornstein-Uhlenbeck process: An extension, European Finance Review 

3 23 (1999). 

[13] J. Masoliver and J. Perello, A correlated stochastic volatility model measuring leverage and other stylized facts, International 
Journal of Theoretical and Applied Finance 5 541 (2002). 

[14] A. Dragulescu and V. Yakovenko, Probability distribution of returns in the Heston model with stochastic volatility, Quan- 
titative Finance 2 443 (2002). 

[15] B. LeBaron, Stochastic volatility as a simple generator of apparent financial power laws and long memory, Quantitative 
Finance 1 621 (2001). 

[16] A. Fisher, L. Calvet, B. B. Mandelbrot, Multifractality of DEM/S rates, Cowles Foundation Discussion Paper 1165. 
[17] L. Calvet, A. Fisher, Multifractality in asset returns: theory and evidence, Forthcoming in Review of Economics and 
Statistics. 

[18] F. Black, Studies of stock market changes, Proceedings 1976 of the American Statistical Association, Business and Eco- 
nomics Statistics Section 177 (1976). 

[19] G. Bekaert, G. Wu, Asymmetric volatility and Risk in Equity markets, The Review of Financial Studies 13 1 (2000). 

[20] J.-P. Bouchaud, A. Matacz, and M. Potters, Leverage effect in financial markets: The retarded volatility model, Physical 
Review Letters 87 228701 (2001). 

[21] J. Perello and J. Masoliver, Random diffusion and leverage effect in financial markets, forthcoming in Physical Review E 
(2003). 

[22] B. Pochart and J.-P. Bouchaud, The skewed multifractal random walk with applications to option smiles, Quantitative 

Finance 2 303 (2002). 
[23] E. Bacry, J. F. Muzy, in preparation. 



19 



[24] E. Bacry, J. Delour, J. F. Muzy, Modelling financial time series using multifractal random walks, Physica A 299 84 (2001). 
[25] S. Micciche, G. Bonanno, F. Lillo, R. N. Mantegna, Volatility in financial markets: stochastic models and empirical results, 
e- print cond-mat/0202527 

[26] J. -P. Bouchaud and M. Potters Theory of Financial Risks (Cambridge University Press, Cambridge, 2000), Theory of 

Financial Risk and Derivative Pricing (Cambridge University Press, Cambridge, 2003, forthcoming). 
[27] R. Engle and A. Patton, What good is volatility model?, Quantitative Finance 1 237 (2001). 

[28] J. Perello, J. M. Porra, M. Montero, and J. Masoliver, Black-Scholes option pricing within It and Stratonovich conventions, 
Physica A 278 260 (2000). 

[29] T. Lux, The stable Paretian hypothesis and the frequency of large returns: an examination of major German stocks, Applied 

Financial Economics 6 463 (1996). 
[30] F. Longin, The Asymptotic Distribution of Extreme Stock Market Returns, Journal of Business, 69 383 (1996) 
[31] V. Plerou, P. Gopikrishnan, L. A. Amaral, M. Meyer, H.E. Stanley, Scaling of the distribution of price fluctuations of 

individual companies, Physical Review E 60 6519 (1999). 
[32] J. Masoliver, M. Montero, J. M. Porra, A dynamical model describing stock market price distributions, Physica A 283, 

559-567 (2000). 

[33] J. Masoliver, M. Montero, G. H. Weiss, Continuous-time random-walk model for financial distributions, forthcoming in 
Physical Review E (2003). 

[34] E. A. Novikov, Functionals and the random-force method in turbulence theory, Soviet Physics JETP 20, 1290 (1965). 
[35] P. Hanggi in Stochastic Processes Applied to Physics, L. Pesquera and M. A. Rodriguez Eds., 69-95, (World Scientific 
Press, Singapore, 1985). 



